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Abstract 

Based on 88 million T(4S') BB decays collected by the BABMi experiment at the PEP- 
II asymmetric-energy B factory at SLAC, we report preliminary results of four analyses which 
investigate semileptonic charmless B decays, B Xu^v. Deeper understanding of all aspects 
of these decays will improve the determination of the Cabibbo-Kobayashi-Maskawa matrix ele- 
ment \Vub\- In events in which one B meson decay to a hadronic final state is fully reconstructed, 
the semileptonic decay of the second B meson is identified by the detection of a charged lep- 
ton. By measuring the spectrum of the invariant mass of the hadronic system {Mx)-, we de- 
rive the branching fraction B{B XJu) = (2.53 ± 0.29(stat.) ± 0.26(sys.)lo:4?(theo.)) x lO'^. 
The two-dimensional distribution of Mx and q'^, the squared lepton- neutrino invariant mass, 
is used to derive the partial branching fraction for Mx < 1.7GeV/c^,g^ > 8GeV^/c^ to be 
AB{B Xulv) = (0.88 ± 0.14(stat.) ± 0.13(sys.) ± 0.02(theo.)) x IQ-^. From these two 
measurements we can extract \Vub\ = (4.77 ± 0.28(stat.) ± 0.28(sys.)lo;45(theo.)) x 10~^ and 
\Vub\ = (4.92 ± 0.39(stat.) ± 0.36(sys.) ± 0.46(theo.)) x lO'^, respectively We use the same 
sample to extract the true Mx distribution for B Xu^u events, with the goal of comparing 
it with theoretical models. We also identify several exclusive charmless semileptonic B decays, 
and measure the branching fractions B{B irlv) = (1.08 it 0.28(stat.) it 0.16(sys.)) x 10"'^ and 
B{B p£u) = (2.57 ± 0.52(stat.) ± 0.59(sys.)) x 10~^ using isospin and quark model constraints. 
We also set limits on B{B- r/£z? ), B{B- ri'lD ), B{B- a^^v )i3(a[] ^ r/7r°), and 
B(B^ a^lD )B{a+ r]Tr+). 

Submitted to the 32°^^ International Conference on High-Energy Physics, ICHEP 04, 
16 August— 22 August 2004, Beijing, China 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



The BABAR Collaboration, 

B. Aubert, R. Barate, D. Boutigny, F. Coudcrc, J.-M. Gaillard, A. Hicheur, Y. Karyotakis, J. P. Lees, 

V. Tisscrand, A. Zghichc 

Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 

A. Palano, A. Pompili 
Universitd di Bari, Dipartimento di Fisica and IN FN, 1-70126 Bari, Italy 

J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigen, I. Ofte, B. Stugu 
University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 

G. S. Abrams, A. W. Borgland, A. B. Breon, D. N. Brown, J. Button-Shafer, R. N. Cahn, E. Charles, 

C. T. Day, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadel, J. Kadyk, L. T. Kerth, 
Yu. G. Kolomensky, G. Kukartsev, G. Lynch, L. M. Mir, P. J. Oddone, T. J. Orimoto, M. Pripstein, 

N. A. Roe, M. T. Ronan, V. G. Shelkov, W. A. Wenzel 

Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720, USA 

M. Barrett, K. E. Ford, T. J. Harrison, A. J. Hart, C. M. Hawkes, S. E. Morgan, A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

M. Fritsch, K. Goetzen, T. Held, H. Koch, B. Lewandowski, M. Pehzaeus, M. Steinke 
Ruhr Universitdt Bochum, Institut fiir Experimentalphysik 1, D-44'if80 Bochum, Germany 

J. T. Boyd, N. Chevalier, W. N. Cottingham, M. P. Kelly T. E. Latham, F. F. Wilson 
University of Bristol, Bristol BS8 ITL, United Kingdom 

T. Cuhadar-Donszelmann, C. Hearty, N. S. Knecht, T. S. Mattison, J. A. McKenna, D. Thiessen 
University of British Columbia, Vancouver, BC, Canada VST IZl 

A. Khan, P. Kyberd, L. Teodorescu 
Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 

A. E. Blinov, V. E. Blinov, V. P. Druzhinin, V. B. Golubev, V. N. Ivanchenko, E. A. Kravchenko, 
A. P. Onuchin, S. L Serednyakov, Yu. L Skovpen, E. P. Solodov, A. N. Yushkov 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. Best, M. Bruinsma, M. Chao, L Eschrich, D. Kirkby, A. J. Lankford, M. Mandelkern, R. K. Mommsen, 

W. Roethel, D. P. Stoker 

University of California at Irvine, Irvine, CA 92697, USA 

C. Buchanan, B. L. Hartfiel 
University of California at Los Angeles, Los Angeles, CA 90024, USA 

S. D. Foulkes, J. W. Gary, B. C. Shen, K. Wang 
University of California at Riverside, Riverside, CA 92521, USA 



2 



D. del Re, H. K. Hadavand, E. J. Hill, D. B. MacFarlane, H. P. Paar, Sh. Rahatlou, V. Sharma 
University of California at San Diego, La Jolla, CA 92093, USA 

J. W. Berryhill, C. Campagnari, B. Dalimes, O. Long, A. Lu, M. A. Mazur, J. D. Richman, W. Verkerke 
University of California at Santa Barbara, Santa Barbara, CA 93106, USA 

T. W. Beck, A. M. Eisner, C. A. Heusch, J. Kroseberg, W. S. Lockman, G. Nesom, T. Schalk, 
B. A. Schumm, A. Seiden, R Spradlin, D. C. Williams, M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, CA 95064, USA 

J. Albert, E. Chen, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. Hitlin, I. Narsky, T. Piatenko, 

F. C. Porter, A. Ryd, A. Samuel, S. Yang 

California Institute of Technology, Pasadena, CA 91125, USA 

S. Jayatilleke, G. Mancinelli, B. T. Meadows, M. D. Sokoloff 
University of Cincinnati, Cincinnati, OH 45221, USA 

T. Abe, F. Blanc, P. Bloom, S. Chen, W. T. Ford, U. Nauenberg, A. Olivas, P. Rankin, J. G. Smith, 

J. Zhang, L. Zhang 

University of Colorado, Boulder, CO 80309, USA 

A. Chen, J. L. Harton, A. Soffer, W. H. Toki, R. J. Wilson, Q. Zeng 
Colorado State University, Fort Collins, CO 80523, USA 

D. Altenburg, T. Brandt, J. Brose, M. Dickopp, E. Feltresi, A. Hauke, H. M. Lacker, R. Miiller-Pfefferkorn, 
R. Nogowski, S. Otto, A. Petzold, J. Schubert, K. R. Schubert, R. Schwierz, B. Spaan, J. E. Sundermann, 

K. Tackmann 

Technische Universitdt Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, G. R. Bonneaud, F. Brochard, P. Grenier, S. Schrcnk, Ch. Thicbaux, G. Vasileiadis, M. Verderi 

Ecole Poly technique, LLR, F-91128 Palaiseau, France 

D. J. Bard, P. J. Clark, D. Lavin, F. Muheim, S. Playfer, Y. Xie 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti, V. Azzolini, D. Bettoni, C. Bozzi, R. Calabroso, G. Cibinetto, E. Luppi, M. Negrini, 

L. Picmontcsc, A. Sarti 

Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

E. TrcadwcU 

Florida A&M University, Tallahassee, FL 32307, USA 

F. Anulli, R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, P. Patteri, I. M. Peruzzi, 

M. Piccolo, A. Zallo 

Laboratori Nazionali di Frascati dellTNFN, 1-00044 Frascati, Italy 

A. Buzzo, R. Capra, R. Contri, G. Crosetti, M. Lo Vetere, M. Macri, M. R. Monge, S. Passaggio, 

C. Patrignani, E. Robutti, A. Santroni, S. Tosi 

Universitd di Genova, Dipartimento di Fisica and INFN, I-I6I46 Genova, Italy 



3 



S. Bailey, G. Brandenburg, K. S. Chaisanguanthum, M. Morii, E. Won 
Harvard University, Cambridge, MA 02138, USA 

R. S. Dubitzky, U. Langenegger 
Universitdt Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

W. Bhimji, D. A. Bowerman, P. D. Dauncey, U. Egede, J. R. Gaillard, G. W. Morton, J. A. Nash, 

M. B. Nikolich, G. R Taylor 

Imperial College London, London, SW7 2AZ, United Kingdom 

M. J. Charles, G. J. Grenier, U. MaUik 
University of Iowa, Iowa City, lA 52242, USA 

J. Cochran, H. B. Crawley, J. Lamsa, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin, J. Yi 
Iowa State University, Ames, I A 50011-3160, USA 

M. Biasini, R. Covarelli, M. Pioppi 
Universitd di Perugia, Dipartimento di Fisica and INFN, 1-06100 Perugia, Italy 

M. Davier, X. Giroux, G. Grosdidier, A. Hocker, S. Laplace, F. Le Diberder, V. Lcpeltier, A. M. Lutz, 
T. C. Petersen, S. Plaszczynski, M. H. Schune, L. Tantot, G. Wormser 

Laboratoire de V Accelerateur Lineaire, F-91898 Orsay, France 

C. H. Cheng, D. J. Lange, M. C. Simani, D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, CA §4550, USA 

A. J. Bevan, C. A. Chavez, J. P. Coleman, I. J. Forster, J. R. Fry, E. Gabatlnder, R. Garnet, 
D. E. Hutchcroft, R. J. Parry, D. J. Payne, R. J. Sloane, C. Touramanis 

University of Liverpool, Liverpool L69 72E, United Kingdom 

J. J. Back,i C. M. Cormack, P. F. Harrison,^ F. Di Lodovico, G. B. Mohanty^ 

Queen Mary, University of London, El 4NS, United Kingdom 

C. L. Brown, G. Cowan, R. L. Flack, H. U. Flaecher, M. G. Green, P. S. Jackson, T. R. McMahon, 

S. Ricciardi, F. Salvatore, M. A. Winter 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, 

United Kingdom 

D. Brown, C. L. Davis 
University of Louisville, Louisville, KY 40292, USA 

J. AlUson, N. R. Barlow, R. J. Barlow, P. A. Hart, M. C. Hodgkinson, G. D. Lafferty, A. J. Lyon, 

J. C. Wilhams 

University of Manchester, Manchester M13 9PL, United Kingdom 

A. Farbin, W. D. Hulsbergen, A. Jawahery, D. Kovalskyi, C. K. Lae, V. Lillard, D. A. Roberts 
University of Maryland, College Park, MD 20742, USA 
^Now at Department of Physics, University of Warwick, Coventry, United Kingdom 



4 



G. Blaylock, C. Dallapiccola, K. T. Flood, S. S. Hertzbach, R. Kofler, V. B. Koptchev, T. B. Moore, 

S. Saremi, H. Staengle, S. Willocq 

University of Massachusetts, Amherst, MA 01003, USA 

R. Cowan, G. ScioUa, S. J. Sekula, F. Taylor, R. K. Yamamoto 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, MA 02139, USA 

D. J. J. Mangeol, R M. Patel, S. H. Robertson 
McGill University, Montreal, QC, Canada H3A 2T8 

A. Lazzaro, V. Lombardo, F. Palombo 
Universitd di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Kroeger, J. Reidy, D. A. Sanders, D. J. Summers, 

H. W. Zhao 

University of Mississippi, University, MS 38677, USA 

S. Brunet, D. Cote, P. Taras 
Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, QC, Canada H3C 3J7 

H. Nicholson 

Mount Holyoke College, South Hadley, MA 01075, USA 

N. Cavallo,^ F. Fabozzi,^ C. Gatto, L. Lista, D. Monorchio, P. Paolucci, D. Piccolo, C. Sciacca 
Universitd di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

M. Baak, H. Bulten, G. Raven, H. L. Snoek, L. Wilden 

NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, 

The Netherlands 

C. P. Jcssop, J. M. LoSccco 
University of Notre Dame, Notre Dame, IN 46556, USA 

T. Allmendinger, K. K. Gan, K. Honscheid, D. Hufnagel, H. Kagan, R. Kass, T. PuUiam, A. M. Rahimi, 

R. Ter-Antonyan, Q. K. Wong 

Ohio State University, Columbus, OH 43210, USA 

J. Brau, R. Frey, O. Igonkina, C. T. Potter, N. B. Sinev, D. Strom, E. Torrence 
University of Oregon, Eugene, OR 97403, USA 

F. Colecchia, A. Dorigo, F. Galeazzi, M. Margoni, M. Morandin, M. Posocco, M. Rotondo, F. Simonetto, 

R. Stroih, G. Tiozzo, C. Voci 

Universitd di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauveau, P. David, Ch. de la Vaissiere, L. Del Buono, O. Hamon, 
M. J. J. John, Ph. Leruste, J. Malcles, J. Ocariz, M. Pivk, L. Roos, S. T'Jampens, G. Therin 

Universites Paris VI et VII, Laboratoire de Physique Nucleaire et de Hautes Energies, F-75252 Paris, 

France 

■^Also with Universita della Basilicata, Potenza, Italy 



5 



p. F. Manfrcdi, V. Re 
Universitd di Pavia, Dipartimento di Elettronica and INFN, 1-27100 Pavia, Italy 



P. K. Bchera, L. Gladncy, Q. H. Guo, J. Panetta 
University of Pennsylvania, Philadelphia, PA 19104, USA 

C. Angelini, G. Batignani, S. Bcttarini. M. Bondioli, F. Bucci, G. Caldcrini. M. Carpinclli, F. Forti, 
M. A. Giorgi, A. Lusiani, G. Marchiori, F. Martinez- Vidal,^ M. Morganti, N. Neri, E. Paoloni, M. Rama, 

G. Rizzo, F. Sandrelli, J. Walsh 

Universitd di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56127 Pisa, Italy 

M. Haire, D. Judd, K. Paick, D. E. Wagoner 
Prairie View A&M University, Prairie View, TX 774.4-6, USA 

N. Danielson, P. Elmer, Y. P. Lau, C. Lu, V. Miftakov, J. Olsen, A. J. S. Smith, A. V. Telnov 
Princeton University, Princeton, NJ 08544, USA 

F. Bellini, G. Cavoto,'* R. Faccini, F. Ferrarotto, F. Ferroni, M. Gaspero, L. Li Gioi, M. A. Mazzoni, 
S. Morganti, M. Pierini, G. Piredda, F. Safai Tehrani, C. Voena 

Universitd di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

S. Christ, G. Wagner, R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, N. De Greet, B. Pranek, N. I. Geddes, G. P. Gopal, E. O. Olaiya 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, S. Emery, A. Gaidot, S. F. Ganzhur, P.-F. Giraud, G. Hamel de Monchenault, W. Kozanecki, 
M. Legendre, G. W. London, B. Mayer, G. Schott, G. Vasseur, Ch. Yeche, M. Zito 

DSM/Dapnia, CEA/Saclay, F-91191 Gif-sur-Yvette, France 

M. V. Purohit, A. W. Weidemann, J. R. Wilson, F. X. Yumiceva 
University of South Carolina, Columbia, SC 29208, USA 

D. Aston, R. Bartoldus, N. Berger, A. M. Boyarski, O. L. Buchmueller, R. Glaus, M. R. Convery, 
M. Cristinziani, G. Dc Nardo, D. Dong, J. Dorfan, D. Dujmic, W. Dimwoodic, E. E. Elscn, S. Fan, 
R. C. Field, T. Glanzman, S. J. Gowdy, T. Hadig, V. Halyo, C. Hast, T. Hryn'ova, W. R. Innes, 
M. H. Kelsey, P. Kim, M. L. Kocian, D. W. G. S. Leith, J. Libby, S. Luitz, V. Luth, H. L. Lynch, 
H. Marsiskc, R. Mcssncr, D. R. Muller, C. P. O'Grady, V. E. Ozcan, A. Pcrazzo, M. Perl, S. Pctrak, 
B. N. Ratcliff, A. Roodman, A. A. Salnikov, R. H. Schindler, J. Schwiening, G. Simi, A. Snyder, A. Soha, 
J. Stelzer, D. Su, M. K. Sullivan, J. Va'vra, S. R. Wagner, M. Weaver, A. J. R. Weinstein, 
W. J. Wisniewski, M. Wittgen, D. H. Wright, A. K. Yarritu, C. C. Young 

Stanford Linear Accelerator Center, Stanford, CA 94309, USA 

P. R. Burchat, A. J. Edwards, T. \. Meyer, B. A. Petersen, C. Roat 

Stanford University, Stanford, CA 94305-4060, USA 

''Also with IFIC, Institute dc Fi'sica Corpuscular, CSIC-Universidad de Valencia, Valencia, Spain 
^Also with Princeton University, Princeton, USA 



6 



S. Ahmed, M. S. Alam, J. A. Ernst, M. A. Saccd, M. Saleem, F. R. Wappler 

State University of New York, Albany, NY 12222, USA 

W. Bugg, M. Krishnamurthy, S. M. Spanier 
University of Tennessee, Knoxville, TN 37996, USA 

R. Eckmann, H. Kim, J. L. Ritchie, A. Satpathy, R. F. Schwitters 
University of Texas at Austin, Austin, TX 18112, USA 

J. M. Izen, I. Kitayama, X. C. Lou, S. Ye 
University of Texas at Dallas, Richardson, TX 15083, USA 

F. Bianchi, M. Bona, F. Gallo, D. Gamba 
Universita di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

L. Bosisio, C. Cartaro, F. Cossutti, G. Delia Ricca, S. Dittongo, S. Grancagnolo, L. Lanceri, P. Poropat,^ 

L. Vitale, G. Vuagnin 

Universita di Trieste, Dipartimento di Fisica and INFN, 1-34121 Trieste, Italy 

R. S. Panvini 
Vanderbilt University, Nashville, TN 31235, USA 

Sw. Banerjee, C. M. Brown, D. Fortin, P. D. Jackson, R. Kowalewski, J. M. Roney, R. J. Sobie 
University of Victoria, Victoria, BC, Canada V8W 3P6 

H. R. Band. B. Cheng, S. Dasu, M. Datta, A. M. Eichcnbaum, M. Graham, J. J. Hollar, J. R. Johnson, 
P. E. Kutter, H. Li, R. Liu, A. Mihalyi, A. K. Mohapatra. Y. Pan, R. Prepost, P. Tan, J. H. von 

Wimmersperg-Toellcr, J. Wu, S. L. Wu, Z. Yu 

University of Wisconsin, Madison, WI 53106, USA 

M. G. Greene, H. Neal 
Yale University, New Haven, CT 06511, USA 



® Deceased 



7 



1 Introduction 



The principal physics goal of the BABAR experiment is to establish CP violation in B mesons and to 
test whether the observed effects are consistent with the predictions of the Standard Model (SM). 
CP violating effects result in the SM from an irreducible phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix which describes the couplings of the charged weak current to quarks. An improved 
determination of the absolute value of the matrix element \Vub\, the coupling strength of the b quark 
to the u quark, will contribute critically to tests of the consistency of the angles of the unitarity 
triangle of the CKM matrix. 

The extraction of |V^b| is a challenge both for theory and experiment. Experimentally, the main 
problem is the separation of 6 — > uif decays from the more abundant b ciu decays. Selection 
criteria applied to achieve this separation generally make the theoretical extrapolation to the full 
decay rate more difficult. Theoretically, inclusive semileptonic rates can be calculated reliably at 
the parton level. However, the dynamics of B meson decays depend on the b quark mass and 
its motion inside the meson. Calculations of the decay rate rely on operator product expansions 
(OPE) in inverse powers of the b quark mass. These depend on the choice of renormalization 
scale and include non-perturbative contributions, resummed into the so-called shape function (SF), 
which introduce model uncertainties. Exclusive branching fractions can also be related to \Vub\, 
although with large errors. In addition, the as yet poorly known dynamics of the decays introduce 
large uncertainties in the determination of the efficiencies of the selection criteria. To reduce these 
uncertainties we apply as loose selection criteria as possible to exclusive B —>■ Xulv decays^. 

At the T(4S') resonance, events with a reconstructed B decay to hadronic final states are an 
optimal environment for the study of semileptonic decays of the second B meson. The relatively 
small backgrounds allow for loose selection criteria and reduce the uncertainty in the extrapolation 
to the full decay rate. Moreover, the fact that the B momentum is known allows us to isolate 
the signal in several regions of phase space and perform different measurements with relatively 
uncorrelated theoretical and experimental uncertainties. The kinematic variables considered in 
this paper are Mx, the invariant mass of the hadrons X, and q^, the squared invariant mass of the 
two leptons. The main background from B Xc^v decays is located at high values of Mx and low 
values of (f, while the B — XuiU signal events extends to low Mx and high values (see Fig 1). 

In this paper we present several preliminary results that extend the |V^b| analysis published 
by BABAR, in Rcf. [1]. These new studies are based on the same analysis strategy and data set. 
They are motivated by the recent discussion on the theoretical uncertainties to be assigned to the 
measurement of \Vub\, see for instance [2, 3, 4] and by the possibility of studying exclusive charmless 
decays with high purity. This paper presents four separate results: 

• A \Vub\ measurement is performed utilizing the Mx spectrum, as in Ref. [1], following the 
approach from De Fazio and Neubert [5] referred to as "DFN" in the following. An updat cKl 
estimate of the SF parameters is used for the extrapolation to the full phase space. This 
analysis will be referred to as the "M^" analysis. 

• A determination of the true Mx spectrum with detector effects unfolded ("unfolding" anal- 
ysis). This spectrum allows for a direct comparison with theoretical models and with more 
statistics will be able to constrain the SF parameters. We determine the cumulative Mx 
distribution, as well as the first and second mass moment of the Mx distribution. 

^Unless otherwise specified, charge conjugation is always imphed throughout this paper. 
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Mx(GeV/c^) Mj^CGeV/c^) 



Figure 1: Distribution on MonteCarlo simulated events of the squared invariant mass of the leptons 
{(f') and the invariant mass of the hadronic recoil system {Mx) in semileptonic b — > ulv (left) and 
h — > civ decays (right). The model utilized is the non-resonant model described in Sec. 2.1. The 
full (dashed) line indicates the phase-space region selected by the Mx{Mx-q^) analysis. 

• A \Vub\ measurement is performed by utilizing the two-dimensional distribution of Mx and 

and by following the approach of Bauer, Ligeti and Luke [6], referred to as "BLL" in 
the following. This analysis, which will be referred to as the "Mx-g^" analysis, requires the 
determination of the partial branching fraction in limited regions of the phase space. 

• A measurement of branching fractions for exclusive charmless semileptonic B decays. In 
this case the semileptonic decay is identified by a charged lepton, and the hadronic state is 
exclusively reconstructed. We analyze the following decay modes: ir'^iu , B~ — > ir^iiy , 
5° ^ p+ii> , B- p^lv , B- -> ujiu , B- rjii? , B' iv , B' a^lv and W 

A review of the previous measurements of \Vub\ is given in Ref. [4]. The measurements reported 
here are either based on novel techniques or have smaller uncertainties than the existing ones. 

This paper is organized as follows: section 2 describes the detector, the data sample and the 
Monte Carlo simulation including a description of the theoretical model on which our efficiency 
calculations are based. Section 3 describes the aspects common to all analyses presented here. 
Section 4 presents a new interpretation of the published Mx analysis that takes the latest theoretical 
developments into account. Section 5 discusses the results of the Mx spectrum unfolding (on the 
whole range). Section 6 shows the details and the results of the Mx-q^ analysis, while Section 7 
shows the results for an alternative set of SF parameters. Finally in Section 8 we present the 
measurement of exclusive B — > Xuiv branching fractions. 
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2 Data Sample and Simulation 



The data used in this paper were recorded with the BABAR detector [7] at the PEP-II colhder 
in the period October 1999-September 2002. The total integrated luminosity of the data set is 
81.9 fb"^ collected on the T{AS) resonance. The corresponding number of produced BB pairs is 
88 million. We use Monte Carlo (MC) simulations of the BABAR detector based on GEANT [8] to 
optimize selection criteria and to determine signal efficiencies and background distributions. 

2.1 Simulation oi B ^ Xu^v decays 

Charmless semileptonic B X^iv decays are simulated as a combination of both three-body 
decays to narrow resonances, X^ = Tr,rj, p,uj,r]' , and decays to non-resonant hadronic final states 




Figure 2: Ax^ = 1 contours for the fit to 6 —> cii' moments in Ref. [9] (green) and the b ^ s'y 
photon energy spectra in Belle [12] (red) and CLEO [4] (blue). Dots represent the best points. 

The simulation of the inclusive charmless semileptonic B decays into hadronic states with masses 
larger than Im-^ is based on a prescription by De Fazio and Neubert [5] (DFN), which calculates 
the triple differential decay rate, d^T / dq^ dE^dsH {sh = "t-x)' ^{(^s) corrections. The 

motion of the h quark inside the B meson is incorporated in the DFN formalism by convolving 
the parton-level triple differential decay rate with a non-perturbative shape function (SF). The SF 
describes the distribution of the momentum /c+ of the b quark inside the B meson. The two free 
parameters of the SF are A^^ and Af^. The first relates the B meson mass, m^, to the b quark 
mass, mf^ = tub — A^^ ^ and — Af''^ is the average momentum squared of the b quark in the B 
meson. The SF parameterization used in the generator is of the form 

F{k+) = N{l-xYe^^+''>, (1) 
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Figure 3: Model for the simulation of 6 ^ uiv' decays for neutral (top) and charged (bottom) B 
mesons: resonant decays (purple - peaky) are combined with weighted non-resonant ("inclusive") 
simulated events (red - smooth) to form the "hybrid" model (black). The inclusive sample is also 
shown before applying the weights (blue). 

where x = < 1 and a = —3{A^^)'^/Xf^ — 1. The original DFN paper [5] suggested that the 
SF parameters could be related to the operator product expansion parameters by — A^^ = 
and Af^ = — /i^. Since there is currently no consensus on these relationships (see for instance 
Refs. [2, 3]) we choose to use the values of A^^ and Af ^extracted from the b ^ s'y spectrum by 
CLEO [4]. Fitting the spectrum to simulated samples with different SF parameters allows us to 
create a Ax^ = 1 contour, shown in Figure 2, which we use to estimate theoretical uncertainties. 
The point with minimum x 

corresponds to A^^ = 0.545 GeV/c^ and Af^ = -0.342 GeV Vc'^- 
In the simulation the hadron system is produced with a non-resonant and continuous invari- 
ant mass spectrum according to the DFN model. A reweighting of the Fermi motion distribution 
is used to obtain distributions for different values of A^^ and Af^. Finally, the fragmentation of 
the Xu system into final state hadrons is performed by JETSET [10]. The exclusive charmless 
semileptonic decays are simulated using the ISGW2 model [11]. The resonant and non-resonant 
components are combined such that the total branching fraction is consistent with the measured 
value [1] and that the integrated spectrum agrees with the prediction of Ref. [5]. The resulting nix 
distributions for charged and neutral B mesons are shown in Fig. 3. All branching fractions and 
theory parameters involved in this reweighting are varied within their errors in the evaluation of 
the associated uncertainty. 

While we were writing this paper we became aware of a preliminary interpretation of the 6 — > 57 
energy spectrum measured by Belle [12]. The resulting Ax^ = 1 contour is shown in Fig. 2: it is 
consistent with the CLEO one and has smaller uncertainties. A visible shift of the central value of 
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the measured quantities is expected though because the best fit result from Behe is outside the la 
CLEO contour. Since there was not enough time to combine the two results, from Belle and CLEO, 
we prefer to report the results with the CLEO ellipse as primary results and quote the results with 
the Belle ellipse as a reference for future developments. 

2.2 Theoretical Acceptance Estimates 

The Mx-q^ analysis results in the partial branching fraction for charmless semileptonic decays in 
selected phase space regions. To translate it into a measurement of the total branching fraction, 
and therefore \Vub\, we need the fraction of events inside the measurement region (referred to as 
"acceptance" in the rest of the paper) as an external input. Wc choose the results of Bauer, 
Ligeti and Luke [6] for the acceptance corrections. These authors perform an operator-product- 
expansion-based calculation which includes C(a^) and 0(1/ ml) corrections in a region of phase 
space where the non-perturbative effects due to the shape function are small. Their results show 
that the theoretical uncertainties due to the extrapolation to the full phase space can be significantly 
reduced with respect to a measurement based on a single Mx cut by moving the Mx cut to the 
highest practical value allowed by the charm semileptonic background, and by reducing the cut 
as low as possible. As a cross-check, acceptances computed with the DFN model presented in the 
previous section are also used in the Mx-q^ analysis. 

2.3 Simulation of Background Processes 

In order to estimate the shape of the background distributions we make use of simulations of 
e^e~ — > T(4S') BB with the B mesons decaying inclusively. The most relevant backgrounds 
are due to 5 — > Xciv events. The simulation of these processes uses an HQET parametrization of 
form factors for B D*eu [13], and models for B D^*'>tt£U [14], and B Deu,D**iU [11]. 

3 Common Analysis Aspects 

The event selection and reconstruction and the measurements of branching fractions follow closely 
the strategy described in Ref. [1] and represents the common base for the analyses presented here. 
However, each analysis may differ in some details. This is particularly true for the measurement 
of exclusive branching fractions. This section describes aspects in common to all studies. Details 
specific to each analysis are provided in the following sections. 

3.1 Event Reconstruction and Selection 

In this paper we study the recoil of fully reconstructed B in hadronic decay modes (Breco), which 
is a moderately pure sample with known flavor and four-momentum. Wc select Breco decays of the 
type B DY, where D refers to a charm meson, and Y represents a collection of hadrons with 
a total charge of ±1, composed of hitt* -|- na-fC^ + nsKg + n47r'^, where ni + n2 < 6, ns < 3, and 
724 < 3. Using D~ and D*~ {D'^ and D*^) as seeds for B^ {B~^) decays, we reconstruct about 1000 
different decay chains. Overall, we correctly reconstruct one B candidate in 0.3% (0.5%) of the 
B^B^ (B^B^) events. The kinematic consistency of a Breco candidate with a B meson decay is 

checked using two variables, the beam-energy-substituted mass mss = \j ^/'^ ~ Pb energy 
difference, = Eb — \/i/2. Here refers to the total energy in the T(4S') center of mass frame. 
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and pb and Eb denote the momentum and energy of the -Breco candidate in the same frame. For 
signal events the rriES distribution peaJiS at the B meson mass, while AE is consistent with zero. 
We require AE = within approximately three standard deviations. 

A semileptonic decay of the other B meson {^B^^^^qh ) is identified by the presence of a charged 
lepton with momentum in the B^gcoil rest frame (p*) higher than IGeV/c . In addition, the 
detection of missing energy and momentum in the event is taken as evidence for the presence of a 
neutrino. The hadronic system X is reconstructed from charged tracks and energy depositions in 
the calorimeter that are not associated with the Breco candidate or the identified lepton. Care is 
taken to eliminate fake charged tracks, as well as low-energy beam-generated photons and energy 
depositions in the calorimeter from charged and neutral hadrons. The neutrino four-momentum pi, 
is estimated from the missing momentum four- vector Pmiss = Pr{4S) ~ P Breco ~ Px — Pt, where all 
momenta are measured in the laboratory frame, and Pr(4S) refers to the T(4S') meson. 

Undetected particles and measurement uncertainties affect the determination of the four- 
momenta of the X system and neutrino, and lead to a large leakage B ^ Xc£i> background 
from the high Mx into the low Mx region. To improve the resolution on these four-momenta this 
analysis exploits the well-known kinematics of the e'^e"^ T(4S) — ^ BB process and performs a 
two-constraint kinematic fit to the whole event. 

In the sample of reconstructed B decays two backgrounds need to be considered: the com- 
binatorial background from BB and continuum events, due to random association of tracks and 
neutral clusters, which does not peak in mss, and the BB background whose m-ES distribution 
has the same shape as the signal. After applying all selection criteria, the remaining combinatorial 
background is subtracted by performing an unbinned likelihood fit to the mss distribution. In 
this fit, the combinatorial background originating from e'^e~ qq {q = u, d, s, c) continuum and 
BB events is described by an empirical threshold function [15], and the signal is described by a 
modified Gaussian [16] peaked at the B meson mass. 

To reject the background in the sample of semileptonic decays we require exactly one charged 
lepton with p* >lGeV/c, a total event charge of zero, and a missing mass consistent with zero 
(''^miss < 0.5GeV^/c^). These criteria partly suppress the dominant B — *■ Xctu decays, many 
of which contain an additional neutrino or an undetected meson. We explicitly veto the 
B^ — > D*~^£~V background by searching candidates for such a decay with a partial reconstruc- 
tion technique, that is only identifying the 7r+ from the D*+ — D'^nf decay and the lepton: since 
the momentum of the is almost coUinear with the D*'^ momentum in the laboratory frame, we 
can approximate the energy of the D*+ as E£)*+ ~ m£)*+ ■ E-j^^/lAbMeV/c^ and estimate the neu- 
trino mass as ml = {ps —Pd*+ ~ PiY ■ Events with > — 3GeV^/c^ are likely to be background 
events and are rejected. Finally, we veto events with charged or neutral kaons in the recoil B to 
reduce the peaJiing background from B Xc(.i> decays. Charged kaons are identified [7] with an 
efficiency varying between 60% at the highest and almost 100% at the lowest momenta. The pion 
misidentification rate is about 2%. The ■k~^'k~ decays are reconstructed with an efficiency of 

80% from pairs of oppositely charged tracks with an invariant mass between 486 and 510 MeV/c^. 

3.2 Measurement of Chcirmless Semileptonic Branching Ratios 

To reduce the systematic uncertainties in the derivation of branching fractions, the observed number 
of signal events, corrected for peaking background and efficiency, is normalized to the total number 
of semileptonic decays B Xiv in the recoil of the Breco candidates. The number of observed 
Breco events which contain a charged lepton with p* > IGeV/c is denoted as N^^"-^. It can be 
related to the true number of semileptonic decays, -/V*r"^ and the remaining peaking background 
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BGsi, estimated with Monte Carlo simulation, by 7V*["^ = {N~ - BG si) / tf ef = N^i/efef . Here 
e^' refers to the efficiency for selecting a lepton from a semileptonic B decay with a momentum 
above pcut in an event with a reconstructed B with efficiency ef. Figure 4 shows the result of the 
rriES fit used to determine N^^""^. 
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Figure 4: Fit to the m-ES distribution for the sample with high momentum lepton. 

If we denote as NJ^^'^^ the number of events fitted in the sample after all requirements, and 
with BGu the peaking background coming from semileptonic decays other than the signal, the true 
number of signal events N^^^ is related to them by 



BGu 



u u jijtrue 
^sel^l H ^^u ! 



(2) 



where the signal efficiency e^^^ accounts for all selection criteria applied on the sample after the 
requirement of a high momentum lepton. 

To measure BGu in the inclusive studies, the peaking background {BGu) is estimated by per- 
forming a fit on the Mx or Mx-Q^ distributions, resulting from mss fits in individual Mx or 
Mx-q^ bins, with the shape of the background estimated from Monte Carlo simulation, and its 
normalization free to vary. An example of such a fit for the Mx analysis is shown in Fig. 5. For 
the exclusive analysis the background normalization is taken from Monte Carlo scaled to the data 
luminosity. 

The ratio between the branching fractions for the signal and B Xlv decays is 



Ru/sl 



S (signal) 

BiB xm 



true 



BGu)/{el 



sell 



Ajtrue 
si 



meas 
si 



BG 



si 



,sl,sl 

^l 



(3) 



The efficiency ratio is expected to be close to, but not equal to unity. Due to the difference in 
multiplicity and the different lepton momentum spectra, we expect the tag efficiency et and lepton 
efficiency ei to be slightly different for the two classes of events. The signal branching fraction is 
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Figure 5: fit to the mx distribution: a) data (points) and fit components, and b) data and 
signal MC after subtraction of the b ciu and the "other" backgrounds. 



then obtained from Ru/si using the B Xc^f branching ratio measured in BABAR, 

B{B XJij) = (10.61 ± 0.16(exp.) ± 0.06(theo.))% [17]. (4) 

which, given the previous BABAR measurement for the charmless semileptonic branching ratio [1], 
corresponds to a total semileptonic branching fraction of B{B Xii?) = (10.83 it 0.19)%. 



3.3 Systematic Uncertainties 

Most of the systematic uncertainties are common to all analyses. Uncertainties related to the 
reconstruction of charged tracks are determined by removing randomly a fraction of tracks corre- 
sponding to the uncertainty in the track finding efficiency (1.3% for tracks with transverse momen- 
tum p± > 0.2GeV/c and 2.5% for tracks with p± < 0.2GeV/c). The systematic error due to the 
reconstruction of neutral particles in the EMC is studied by varying the resolution and efficiency 
to match those found in control samples in data. 

We estimate the systematic error due to particle identification by varying the electron and kaon 
identification efficiencies by ±2% and the muon identification efficiency by ±3%. The misidentifi- 
cation probabilities are varied by 15% for all particles. 

The uncertainty of the Breco combinatorial background subtraction is estimated by changing 
the signal shape function to a Gaussian function instead of the empirical function of Ref. [16]. We 
evaluated the effect of cross-feed between B^ and B~^ decays by repeating the analysis with only 
the B° B° or the B+ B' Monte Carlo. 

The impact of the charm semileptonic branching fractions has been estimated by varying each 
of the exclusive branching fractions within one standard deviation of the current world average [18]. 
Similarly, the branching fractions of charm mesons have been varied, both for exclusive decays and 
for inclusive kaon production. To study the cross-feed among the charmless modes we also varied 
the number of charmless exclusive semileptonic decays by 30% for B — > niu and B plu, by 40% 
for B — > uiu and by 100% for the remaining exclusive charmless semileptonic B decays. 
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In presenting the results the systematic effects are divided into the following categories: un- 
certainties related to detector and reconstruction simulation (tracking reconstruction, neutral re- 
construction, lepton identification, charged Kaon identification), adet', uncertainties related to the 
theoretical model of the non-resonant b uii> decays, cJtheo] uncertainties related to other aspects 
of the h — > uiv simulation ( resonant signal decay branching fractions, ss production, hadronization, 
etc.), auiv', uncertainties due to the limited available Monte Carlo statistics {(TMCstat)\ uncertain- 
ties related to background simulation and subtraction {B —>■ D^*'**^lv and D branching fractions, 
fit to the Mx distribution), a^kg', uncertainties related to the subtraction of the combinatorial 
background and cross-feed, (Jhreco- 

4 One-dimensional Mx analysis 

We already published a measurement of \Vub\ on B ^ XuC-y events with Mx < 1.55GcV/c^ [1] and 
subtracting backgrounds with a one-dimensional fit to the Mx distribution (see Fig. 5). The values 
for the shape function parameters used were A^^ = 0.48 ± 0.12 GeV and Af^ = —0.30 ±0.11 GeV^ 
and a correlation of -80% [9] (see Fig. 2). Since this publication, there have been developments in the 
theoretical interpretation. Here we report the updated measurement of \Vub\ obtained by varying 
the shape function parameters as stated in Sec. 2.1, the new central value being A^^ = 0.545 GeV/ 
and Af^ = —0.342. The experimental systematic errors considered are the same and correspond 
to the sum in quadrature of adet, (^ulu, (^bkg^ f^MCstat, and aiyreco- The measured total inclusive 
branching ratio is 

i?^/^; = 0.0234 ± 0.0027(stat.) ± 0.0026(sys.)l|]:|]g^^ (theo.) (5) 

which translates into a total inclusive branching fraction 

B{B Xuiv) = (2.53 ± 0.29(stat.) ± 0.28(sys.)lgj?(theo.)) x 10"^ (6) 

From the relationship between the inclusive branching fraction and \Vub\ (obtained from Eq. 2 in 
Ref. [4] having updated the OPE parameters to the measurement of BABMi [17]) 

IK.I =0.00424(^^^^11^ X (1.0±0.048(OPE + m,)). (7) 

Here the average B^ and B^ lifetime is tb = 1.604 ps [19]. The ^^OPE + mb" uncertainty is due to 
perturbative and non perturbative corrections , and to the uncertainty on the 6-quark mass (m^). 
We then obtain 

\Vub\ = (4.77 ± 0.28(stat.) ± 0.28(sys.)t|jj^(theo.) ± 0.23(OPE + mb))10-^ (8) 

consistent within the shape function errors with the published |V^b| = (4.62 it 0.28(stat.) it 
0.28(sys.) ± 0.40(theo.) ± 0.26(OPE + m^)) x lO'^. 

5 Unfolding the Mx Spectrum 

The observed Mx spectrum is one of the main results of the Mx analysis (see Fig. 5). In order to 
convert this into a universal observable, we need to unfold detector and selection effects. To increase 
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sensitivity to the true distribution, the spectrum is re-binned in 310 MeV/ (? bins in order to match 
the resolution. Prom the measured Mx spectrum, represented by a vector, 'x^^^ ^ we extract the 
unfolded spectrum, "x*''". The relationship between the two spectra is A ~x'"''°- = "x [20], where 
the detector response matrix A describes the effects of limited efficiency and finite resolution of the 
measurement. A is estimated on signal MC with the hybrid model described in Sec. 2.1 and with 
A^^ and Af ^ set to te CLEO best fit values. To reduce systematic uncertainties the simulation has 
been refined: to account for the unknown details of the fragmentation of quarks into hadrons, we 
divide the sample into several categories according to the event multiplicity of charged and neutral 
particles and adjust the fraction of events in each category in MC to match the one observed in 
data. 

In practice, ^4 is a non-invertible matrix. Additionally, the uncertainties on the measured 

spectrum have to be taken into account adequately. We use an unfolding method based on a singular 
value decomposition of the detector response matrix and a suitable regularization procedure of the 
unfolded result as described in Ref. [20] . 

The systematic uncertainties considered are those described in Sec. 3. The systematic effects 
on the measured spectrum x"^^^ arc evaluated, each giving rise to a covariance matrix, CJ^^^, 
on the measured spectrum. To propagate the covariance matrices for the measured distribution of 
MxiCj^°^) to covariance matrices on the unfolded spectrum C^^^ , two different approaches have 
been taken depending on whether the response matrix is affected by the systematic under study. 
In the case of a constant response matrix, e.g. for background modelling, the observed spectrum is 
smeared a large number of times according to Cj!^^^, the spectra are unfolded and the variations of 
the result in each bin give rise to C^^^. When the systematics affect the response matrix, e.g. for 
tracking and neutral efficiencies, the whole measurement is redone with the appropriate changes in 
reconstruction and the spread serves for the determination of C^^^. The covariance matrix 
is computed as a sum of the individual covariance matrices arising from the different systematic 
errors, C^'^^ = 

The resulting unfolded spectrum needs to be corrected for a small possible bias arising from 
the regularization procedure. To estimate this bias we use the unfolded spectrum Ic^^^, apply the 
detector response matrix A to it and fluctuate the resulting spectrum using the covariance matrix 
on the measured spectrum. These spectra are then unfolded and the bin-by-bin mean deviation 
from the original Ic"^^ is taken as bias. 

5.1 Unfolding results 

Figure 6 shows the unfolded spectrum Ic^^^ normalized to unit area and its cumulative. The 
errors on the spectrum are given by the square-roots of the diagonal elements of the covariance 
matrix C*""^. Detailed listings of "x""^ and C^^^ are given in Table 7 and 8 in Appendix A. 
The uncertainties on the cumulative distribution take the bin-by-bin correlations of the unfolded 
spectrum into account. In Fig. 6 we also show a comparison of the data with Monte Carlo samples 
simulated with SF parameters corresponding to the CLEO best fit and extremes of the Ax^ = 1 
contour. The sensitivity of these spectra to the SF parameters is comparable with the CLEO 
6 — ^ S7 constraint. 

Finally, we compute the first moment Mi = J^'''° Mxf{Mx)dMx /{J^'''° f{Mx)dMx) and 

second central moment M'2 = J^'^^'^Mx - Mif f{Mx)dMx /{J(^'''° f{Mx)dMx) from the un- 
folded spectrum f(Mx), for different upper bounds Mx,o on the Mx spectrum. Table 1 shows the 
measured moments for two different values of Mx,o- 5GeV/c^, corresponding to no cut at all, and 
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Figure 6: The unfolded spectrum normalized to unit area (left) and its cumulative distribution 
(right) as a function of Mx- The spectrum and the cumulative distribution from MC with the 
CLEO best fit A^^ = 0.545 GeV/c^ and Af^^ = -0.342 is shown in green (middle). The orange 
(highest) and violet (lowest) spectra and cumulative distributions correspond to the two extreme 
points of the CLEO ellipse (see Fig. 2), A^^ = 0.800 and 0.435 GeV/c^ and to Af^ = -1.22 and 
—0.16, respectively. In the left plot black errors are only statistical, while the red ones include 
always systematics. 



1.86GeV/c^. The quoted uncertainties take the bin-by-bin correlations of the unfolded spectrum 
into account. With a higher precision of the measurement it will be interesting to use moments of 
the hadronic mass spectrum to constrain the shape function parameters. 



Table 1: First and second central moments of the unfolded Mx spectrum. The correlation between 
the first and second central moment with the same Mx,o is also reported. Finally the breakdown of 
the error into statistical uncertainty (cTstat ©CMCstat), uncertainties related to detector effects (crdet); 
uncertainties related to signal modeling (crui;^ © dthco) and to background modeling and subtraction 
((Tbkg©<7breco) IS given in GeV/c^ and GeV^/c'^, respectively. The values of TWi are in GeV/c^ while 
the values of are in GeV^/c^. 





Mx,o (GeV/c2) 


M 


a{M) 


Correlation 


fstat © fMCstat 


O'dct 


(Julu © 0"theo 


Tbkg © '^brcco 


Ml 


1.86 


1.355 


0.084 




0.064 


0.022 


0.017 


0.023 




1.86 


0.147 


0.034 


-0.819 


0.027 


0.010 


0.005 


0.011 




5 


1.584 


0.233 




0.166 


0.073 


0.060 


0.101 


K 


5 


0.270 


0.099 


0.796 


0.055 


0.037 


0.023 


0.068 



6 Two-dimensional Mx-q^ analysis 

The Mx analysis is systematically limited by the dependence on the shape function. This can be 
overcome by selecting a phase space region where the shape function effects are small, namely the 
region at large values [6]. Therefore we find a trade-off between the statistical and theoretical 
uncertainties by loosening the Mx cut and applying a one. Moreover, since most of the theo- 
retical uncertainties are due to the extrapolation from a selected kinematic region to the full phase 
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Figure 7: Distributions of in bins of Mx- Points are data, the blue, light gray and yellow 
histograms represent the fitted contributions from h — > u^v events with true Mx< 1.7GeV/c^, q^> 
SGeV^c^ , h — > ulv events not satisfying these requirements, and background events, respectively. 

space, measurements of partial branching fractions in different regions of phase space and their 
extrapolation to the full phase space can serve as tests of the theoretical calculations and models. 

In order to extract the partial charmless semileptonic branching ratio in a given region of the 
Mx-q^ plane AB{B — > Xuiv), we define as signal the events with true values of the kinematic 
variables in the chosen region, treating as background those that migrate from outside this region 
because of the resolution. This means that in applying Eq. 3 we include the b — > ulv events outside 
the signal region in BGu and the quoted efficiencies refer only to events generated in the chosen 
{Mx-Q^) region. These efficiencies are computed on Monte Carlo, and therefore are based on the 
DFN model. However, the associated theoretical uncertainty on the final result is small compared 
to the extrapolation error to the full phase space. We divide the events into two-dimensional 
bins of Mx and g^, we fit the rriEs distribution to extract the yield in each bin, and we perform 
a two-dimensional binned fit of the entire Mx-q^ distribution in order to extract the signal and 
background components. The result of the fit is shown in Fig. 7. 

Fig. 8a and Table 2 show, for Mx < 1.7GeV/c^, the partial branching fraction AB{B Xulv) 
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as a function of the minimum cut. 

We convert the measured A.B{B —>■ Xuii^) into \Vub\ by 



I 1927r3 ABjB^XJu) 

'^""'^V^^^ — ^ — 

where tb = 1.61ps and G is a theoretical parameter calculated in the BLL approach [6]. The first 

factor under the square root is 192tt^ /{TsGj^ml) = 0.00779. The measured iKil 

function of 

the cut is shown in Fig. 8b for the acceptances computed by BLL and by DFN. Note that, since 
the operator product expansion breaks down when going to low q^, the BLL calculation is only 
possible for higher values of g^. 

The error on the acceptance as computed by BLL increases for tighter cuts on g^. For smaller 
values of g^, the shape function effects increase. In the signal region q"^ > 8GeV^/c*^, Mx < 
LTGeV/c^ we obtain: 

AB{B ^ XJiy,Mx < L7GeV/c^g^ > 8GeV7c^) = (0.88±0.14(stat.)±0.13(sys.)±0.02(theo.))xl0"^. 

(10) 

To extract we take G as computed by BLL and rescale it to the 6-quark mass as measured 

by BABAR[17], obtaining G = 0.282 ± 0.053, corresponding to an acceptance €bll = 0.325 ± 0.061. 
Eq. 9 yields 

\Vub\ = (4.92 ± 0.39(stat.) ± 0.36(sys.) ± 0.46(theo.)) x 10^^ (11) 

In the DFN model the calculated acceptance is e = 0.337^Qgy4 and by using Equation 7 we 
obtain \ Vub\ = (4.85 ± 0.39(stat.) ± 0.36(sys.) lQ'29(theo.)) x 10"^, in agreement with the extraction 
based on BLL, as well as with the result form the one-dimensional Mx fit. Figure 8 shows the 
measured values for \ Vub\ as a function of the q'^ cut for Mx < 1.7 GeV/c^, showing good consistency 
between the different cuts and theoretical framework. Checks were done also with a looser {Mx < 
1.86 GeV/c^) and a tighter {Mx < 1.5 GeV/c^) cut on Mx, and they give consistent results. 



Table 2: Partial branching fraction measurements (in 10"^ units) for Mx < 1-7 GeV/c^ and 
q^ > q'^y^f, as a function of g^«f "^^^ different sources of uncertainty (as described in 3.3) are also 
reported. 





AB{B XJu) 




O'det 


^breco 


O'bkg 


C^theo 




<^MCstat 





1.68 


0.22 


0.15 


0.08 


0.12 


-0.045 +0.035 


0.137 


0.08 


2 


1.52 


0.20 


0.16 


0.07 


0.11 


-0.028 -F0.036 


0.110 


0.07 


4 


1.33 


0.18 


0.11 


0.06 


0.10 


-0.040 -F0.026 


0.116 


0.06 


6 


1.10 


0.16 


0.14 


0.05 


0.08 


-0.022 -F0.018 


0.083 


0.05 


8 


0.88 


0.14 


0.09 


0.04 


0.06 


-0.028 +0.009 


0.053 


0.05 


10 


0.55 


0.11 


0.03 


0.02 


0.04 


-0.006 +0.019 


0.027 


0.03 


12 


0.41 


0.09 


0.04 


0.02 


0.03 


-0.010 +0.000 


0.033 


0.03 


14 


0.21 


0.06 


0.01 


0.01 


0.02 


-0.012 +0.012 


0.018 


0.02 
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Figure 8: (a) Measured partial branching ratio for Mx < 1.7 GeV/c^ and > Qcut^ ^ function of 
g^^j. The error bar is the sum in quadrature of statistical, systematical and theoretical uncertainties, 
(b) Measured value of \Vub\ for Mx < 1.7GeV/c^ as a function of the cut applied when using 
acceptances from DFN (open points) and BLL (solid points). The error bars include the statistical, 
systematic and theoretical uncertainties, added in quadrature. 
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Figure 9: The unfolded spectrum (left) and its cumulative distribution (right) as a function of Mx- 
The spectrum and the cumulative distribution from MC with the Belle best fit A^^ = 0.66GeV/c^ 
and = — 0.4GeV^/c^ is shown in green. The orange and violet spectra and cumulative distri- 
butions correspond to the two extreme points in the Belle ellipse (see Fig. 2), A^^ = 0.600 GeV/c^ 
and A^P = 0.748 GeV/c2 and to Af^ = -0.66( GeV/c^)^ and Af^ = -0.28( GeV/c^)^ respectively. 
In the left plot black errors are only statistical, while the red ones always include systematics. 

7 Results based on Belle's estimate of the SF parameters 

We report here the results obtained with the SF parameters as estimated from the 5 — > 57 photon 
energy spectrum measured by Belle (see Sec. 2.1) and we reinterpret our results by using them in 
the Mx and Mx-q^ analyses. The acceptance obtained for the DFN model is lower, and therefore 
the charmless semileptonic branching fraction and \ Vuh\ are higher. The theoretical systematics due 
to the shape function parameters is reduced, due to the significantly better precision. 
For the Mx analysis we get 

Ru/si = (2.81 ± 0.32(stat.) ± O.31(sys.)+[5:|(theo.)) x 10"^ (12) 
which translates into 

\Vub\ = (5.22 ± 0.30(stat.) ± 0.31(sys.)lg;^^(SF) ± 0.25(pert + l/mb3))10-^ (13) 

As far as the unfolding is concerned. Fig. 9 compares the measured spectra with the distributions 
corresponding to the SF parameters measured by Belle. 

The partial branching fraction measurements as a function of the q'^ cut obtained by the Mx- 
g^analysis are reported, for Mx< 1.7 GeV/c^, in Table 3. 

The measurement of the partial branching fraction B{B Xulv) in the region limited by 
Mx < 1.7GeV/c2, > QeY/c'^f is 

A13{B ^ XJu,Mx < 1.7GeV/c^g2 > 8GeV7c^) = (0.90 ± 0.14(stat.) ± 0.14(sys.)l|]:[|2(theo.)) x 10"^. 
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By using G = 0.282 ± 0.053 from BLL, wc get 

\Vub\ = (4.98 lb 0.40(stat.) ± 0.39(syst.) ± 0.47(theo.)) x 10"^. 



The DFN acceptance computed at Mx< 1.7 GcV/c^ and q^> 8 ( GeV/c^)^ with the Belle SF 
parameters is e = 0.300]']q'q28- This gives in the DFN framework \Vuh\ = (5.18 zt QAlgtat =t 

Figure 10 shows the results for \Vub\ as a function of the cut for Mx< 1-7 GeV/c^, for both 
DFN and BLL. The two models arc still consistent within the present accuracies. The stability of 
the result and the agreement between the two methods seems to indicate that OPE is still valid in 
this range. 

8 Measurement of Exclusive Charmless Semileptonic Branching 
Fractions 

8.1 Reconstruction of Exclusive Modes. 

The reconstruction and selection of the events follows closely the one described in Sec. 3. Hadrons 
are reconstructed in the following modes: 

• 7r° candidates are defined as pairs of photons with an energy in the laboratory > 30 MeV. 
In order to further reject combinatorial background we apply a cut on the energy of the more 
energetic photon in the recoiling B meson rest frame at E* > 300 MeV. 

• /9° candidates are reconstructed using pairs of charged tracks with opposite charge assuming 
the pion mass. A cut on the momentum of the pions in the recoiling B meson rest frame is 
applied. The tracks are ordered depending on their momentum and we apply a cut on the 
more energetic track at p* > 350MeV/c and on the other one at p* > 150MeV/c. 

• p+ candidates are reconstructed using one charged track (with opposite charge with respect 
to the most energetic lepton in the recoil) and one 7r°. We apply a cut on the momentum of 
the 7r° momentum at p*o > 150MeV/c. 



Table 3: BELLE ellipse: Partial branching fraction AB{B Xuii?) measurements (in 10~^ units) 
for different cuts. Mx is required to be less than 1.7 GeV/ c^. The different sources of uncertainties 
are also reported. 



Qcut > 


AB{B -> Xuiu) 


O'stat 


<^det 




<^bkg 






CMCstat 





1.740 


0.231 


0.159 


0.078 


0.129 


-0.042 +0.026 


0.158 


0.078 


2 


1.584 


0.205 


0.165 


0.071 


0.117 


-0.037 +0.027 


0.145 


0.068 


4 


1.381 


0.186 


0.113 


0.062 


0.102 


-0.036 +0.026 


0.140 


0.061 


6 


1.135 


0.161 


0.144 


0.051 


0.084 


-0.025 +0.017 


0.105 


0.053 


8 


0.896 


0.143 


0.091 


0.040 


0.066 


-0.017 +0.012 


0.064 


0.047 


10 


0.566 


0.113 


0.026 


0.025 


0.042 


-0.006 +0.013 


0.041 


0.036 


12 


0.406 


0.085 


0.038 


0.018 


0.030 


-0.002 +0.003 


0.034 


0.026 


14 


0.207 


0.059 


0.014 


0.009 


0.015 


-0.007 +0.002 


0.026 


0.019 
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Figure 10: BELLE ellipse: Values for \Yub\ as a function of the q'^ cut for Mx< 1.7 GeV/c^ by 
taking acceptances from Bauer, Ligeti and Luke (points) and De Fazio- Neubert (squares). The 
error is the sum in quadrature of the statistical, systematic and theoretical uncertainties. 



• to candidates are reconstructed in the mode to —>■ 7t^it~tt^ {^{^ ~^ vr^vr^vr'^) = 89.1%) by 
using pairs of charged tracks with opposite charge and one vr''. In order to further reject 
combinatorics we apply a cut rejecting combinations that are at the edge of the Dalitz plot 
of the three pions. 

• T] candidates are reconstructed in three decay modes: jy — > 77 {BF = 39.4%), rj — > tt^tt^tt^ 
{BF = 22.6%) and rj tt^tt^tt^ {BF = 32.5%). 

• r]' candidates are reconstructed in two decay modes: rj' pj (BF = 29.5%) and r/' t/tt+tt" 
(BF = 44.3%). We apply cuts on the mass of the daughters at 0.508 GeV/c^ < < 
0.588 GeV/c^ and 0.625 GeV/c^ < M^o < 0.925 GeV/c^. 

• Oq and candidates are reconstructed in the rjir^ mode, ry's are required to have a recon- 
structed mass within the range 0.518 GeV/c^ < < 0.578 GeV/c^. We assume a Breit- 
Wigner width of 71 ± 7MeY/c^ [21] for both a[j and . 

In addition to these criteria we also select mass windows around the nominal meson masses as 
detailed in Table 4 



8.2 Event based cuts 

The event variables used to reject background are the same as described in Sec. 3, but the cut values 
are different and have been optimized mode by mode. Also, the missing mass squared {m'^^gj is 
computed from the lepton and just the charged tracks and photon candidates that form one of 
the specific hadronic final states we are looking for. This means that if more than one mode is 
reconstructed in a given event, more than one value of m^^gg is calculated. For each decay mode 
the best candidate in an event is chosen on the basis of the built from the reconstructed invariant 
mass of the resonance daughter and rrimiss- addition we have applied discriminating variables: 
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the total photon energy {Eneutral)^ the invariant mass of the lepton and the charged track {rritrk), 
the number of reconstructed 7r° and charged tracks in the recoil. A summary of selection criteria 
per mode is shown in Table 4. 



Table 4: Summary of event based cuts and efficiency for all modes. The calculation of the efficiency 
does not include hadronic B reconstruction efficiency (e^). 



Common cuts 


p* > l.OGeV/c, Niepton = 1, Qb{recoil)Qe > , 

Qtot = 0,Nk = 0, 
no additional charged track in the recoil 






B- ir^'iP 


B^' p+iu 


<^.s(GcV7c4) 
Mass (GeV/c2) 
Specific cuts 


\^iniss{T^^)\ < 0-3 

Eneutval < 0.25 GeV 
\mtrk - 3.1| < 0.02 GeV/ c2 


-0.5 < m^,,,(7r^) < 0.7 
0.11 <m< 0.16 


\mmiss{p'^)\ < 0-4 

0.55 < m < 1.0 


Final eS.{esem) 


0.302 ± 0.005 


0.263 ± 0.007 


0.154 ±0.004 




B- p^lu 


B- ojIv 


B- r]ii? 


<.ss(GeV7c^) 
Mass (GeV/c2) 
Specific cuts 


Kis.(p')l<0.4 
0.6 < m < 1.0 

Af(7r0) < 2 

<i..H<-0-4GeV7c4 


l"^missMI <0.4 

0.74 <m< 0.82 
iV(7r°) < 6 


l^missir])] < 0.5 

0.515 <m< 0.575 

<i..(^°)<-l-5GeV7c4 


Final eS.{eseiei) 


0.214 ± 0.006 


0.109 ±0.004 


0.133 ±0.006 






B- ^ a{]f// 




Mass (GcV/c2) 
Specific cuts 


\miniss{v')\ < 0.5 
0.92 <m< 0.99 


\minissH)\ < 0-5 

0.92 < m < 1.04 
mliissi^) > l-5GeV7c4 
<i..('7)>0.5GeV7c4 


\m'Lssiao)\ < 0.5 
0.92 <m< 1.04 
"iLs.(P+)>l-0GeV7c4 


Final eff.(ese;e/) 


0.066 ± 0.004 


0.049 ± 0.008 


0.074 ± 0.007 



8.3 Systematic Uncertainties 

The systematics are summarized in Table 5. The sources common to the other analyses are de- 
scribed in Sec. 3, but there are also sources of systematic uncertainty which are specific to this 

exclusive analysis. 

We evaluate the impact of different form factor calculations changing the ISGW2 model (our 
default) to light cone sum rule calculations by reweighting events. This error is included in a^ii^. 

The systematic effect due to the non-resonant structure has been evaluated by using a different 
mixture that allows the non-resonant contribution to go down to 2tt masses (the same model used 
in the inclusive analysis [1]). The B~ p^lv channel has a different treatment since the tt+tt" 
non-resonant contamination may not be negligible and it cannot be distinguished from the signal 
component. In order to quantify the non-resonant contribution wc use an approach that is very 
similar to the one described in [22]. From measurements of e+e^ data and r decays and from 
Bose symmetry considerations the non-resonant contribution has to come from the isospin / = 
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component. The isospin relationships therefore predict no non-resonant component in p'^ii/ , 

and a ratio 2:1 between the non resonant contribution due to B~ tt'^tt^Iv and B" tt+tt^^P. 
By measuring the B~ — ^ ir^TT^ii? from our data we obtain 

B{B- TT^-K^lu) < 0.6 X 10-^ (90 % C.L.) (14) 

in the mass region 0.6GeV/c^ < m(7r''7r°) < l.OGeV/c^ . The 68% C.L. upper hmit is used to 
evaluate the systematic uncertainty due to the non-resonant component. In the following we will 
quote B{B~ "/9°" ), assuming no non-resonant contribution since we measured a rate for 
B~ TT^TT^iu which is compatible with zero. On the other hand this systematic uncertainty is 
taken into account in the combination of the B^ — > p'^iv , B~ — > p^lv and B~ — > ujiv results 
where isospin and quark model relations are used (see next section). This component is included 

in Gbkg- 



Table 5: Summary of the systematic uncertainties in the measurement of the Ru/si defined in Eq. 3. 



i?«A/[xlO-1 



5^ 



TT 



Uncertainty on ■R„/^;[xlO 



B- 



1^ 



p+iv 



B- 



B' 



0.86 
0.33 



0.81 
0.25 



3.3 
1.1 



0.92 
0.34 



1.12 
0.49 



<^det 

^breco 

(^bkg 

<^MCstat 



0.015 
0.090 

0.043 
0.029 
0.048 



0.063 
0.050 

0.046 
0.028 
0.065 



0.250 
0.257 

0.475 
0.234 
0.257 



0.026 
0.071 
0.047 
0.059 
0.092 



0.076 
0.074 

0.064 
0.078 
0.146 



Total syst. error 



0.115 



0.117 



0.652 



0.141 



0.207 



B- 



r]£u 



B- 



rfiv 



B- 



(^stat 



0.34 
0.36 



2.4 
1.0 



1.8 
1.0 



0.5 
1.3 



^breco 
<^bkg 

CfMCstat 



0.167 
0.046 
0.038 
0.004 
0.073 



0.199 
0.122 
0.106 
0.034 
0.290 



0.754 
0.106 

0.072 
0.180 
0.293 



0.081 
0.133 
0.047 
0.050 
0.204 



Total syst. error 



0.192 



0.389 



0.839 



0.266 



8.4 Results of the exclusive measurements. 

For the extraction of the ratio of branching ratios we use Equation 3. The only difference is 
in the inclusive semileptonic branching ratios since we use B{B^ — Xiu) = (10.4 it 0.3)% and 
B{B~ XiU) = (11.3 lb 0.3)%, obtained from the inclusive semileptonic branching ratio [17] and 
the lifetime ratio between neutral and charged B mesons [18]. The results are summarized in Table 
6. In Figs. 11-19 the projections of the result on the reconstructed mass and on the m^^^g variable 
are shown for each mode. All selection criteria are applied except for the ones on the variable that 
is plotted. 
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Table 6: Fit results for all modes. 





excl 


BGexd 


^excl 












11.1 ±3.9 


0.9 ±0.5 


0.65 ± 0.03 


15350 ± 200 


0.84 ±0.01 


0.86 ± 0.33(stat) 


B- - 




15.5 ±4.1 


2.3 ± 1.2 


0.56 ± 0.03 


25250 ± 300 


0.87 ±0.01 


0.81 ± 0.25(stat) 


B^- 


p+£i? 


20.1 ± 5.7 


2.4 ± 1.7 


0.31 ±0.02 


13100 ± 200 


0.76 ±0.01 


3.3± l.l(stat) 


B- - 




15.7 ±4.3 


4.0 ± 1.0 


0.45 ± 0.03 


22500 ± 200 


0.79 ±0.01 


0.92 ± 0.34(stat) 


B- - 


u)lu 


9.3 ±3.3 


1.7 ±0.8 


0.21 ±0.02 


25250 ± 300 


0.77 ±0.01 


1.12±0.49(stat) 


B- - 


rjiu 


3.8 ±2.7 


1.3 ±0.6 


0.28 ±0.01 


23050 ± 200 


0.87 ±0.02 


0.34 ± 0.36(stat) 


B- - 


r^tv 


13.9 ±4.2 


4.3 ± 1.2 


0.14 ±0.01 


23050 ± 200 


0.87 ±0.02 


2.4± l.O(stat) 


B- - 


allv 


9.1 ±3.5 


2.5 ±0.9 


0.11 ±0.01 


23050 ± 200 


0.91 ±0.03 


2.4± 1.3(stat) 


B^- 


-> a^iv 


3.0 ±3.5 


1.4 ±0.9 


0.16 ±0.02 


14700 ± 200 


0.87 ±0.03 


0.6± 1.4(stat) 
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Figure 11: B^ — > t^^Iu : Projection of the fit result onto the rn^^g^ variable. 
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Figure 12: B tt^£i/ : Projection of the fit result onto the m(7r'') (left) and m^jg^(right) variables. 
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Figure 13: B° 
variables. 



p'^iv : Projection of the fit result onto the m(p"' 
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Figure 14: B 



ih' : Projection of the fit result onto the m{pi~^TT ) (left) and m^-^^ (right) 



variables. Here "p means that no i? — > vr^vr non-resonant contribution is taken into account. 
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15: B~ ujIu : Projection of the fit result onto the m{uj) (left) and m^^^,^ (right) variables. 
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Figure 16: B — > r]£i' : Projection of the fit result onto the m(r/) (left) and m^j^^ (right) variables. 



Since the branching ratios of the Oq , Qq decays are not known, the results for these channels 
will include the Og — > ryvr'^ and t/tt"'" branching ratios. Note nevertheless that they are known 

to be close to unity because they are dominant decay modes. 

We obtain: 



B{B- 
B{B^ 



B{B° - 


-K+lv ) 


= (0.89 ± 0.34(stat.) ± 0.12(sys.)) x lO""^ 


B{B- - 


TT^ii? ) 


= (0.91 ± 0.28(stat.) ± 0.14(sys.)) x lO"'' 


B{B^ - 


p+£i? ) 


= (3.5 ± l.l(stat.) ± 0.7(sys.)) x 10^*^, 


B{B- ' 


'pO" £i? ) 


= (1.04 ± 0.39(stat.) ± 0.16(sys.)) x 10"^ 


B{B- 


ujIv ) 


= (1.26 ± 0.55(stat.) ± 0.24(sys.)) x 10~^ 


B{B~ 


rjlu ) 


= (0.39±0.41(stat.) ± 0.22(sys.)x)10-^, 


B{B- 


r(lv ) 


= (2.7 ± 1.2(stat.) ± 0.5(sys.)) x 10""^, 


^ a^Q^u )B{al 


^7?7rO) 


= (2.7 ± 1.4(stat.) ± 0.9(sys.)) x 10~^, 


a^lv )B{a^ 




= (0.7 ± 1.6(stat.) ± 0.3(sys.)) x 10"^, 



(15) 



where the errors on i3(-B ^ Xiv') are added in quadrature. 

The rates for exclusive decays involving ir'^ and vr*^ and p'^ , , and w can be constrained by the 
following isospin and quark model relations: 



r{B - 

r{B 



niU) ; 

>peu) 



T{B° 
T{B^ 
T{B- 



TT' iU ) 

p+iu ) 



2T{B' 
2T{B- 

r(B- 



p^ei? 
p^liJ ) 



(16) 



and it is possible to combine the results. Using a lifetime ratio of = 1.086 it 0.017 we obtain: 
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Figure 17: B — > r]'iu : Projection of the fit result onto the m(?y') (left) and m^j^^ (right) variables 
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Figure 18: B QqIv : Projection of the fit result onto the m(og ) (left) and m^-^^ (right) variables. 
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Figure 19: S° 
variables. 



: Projection of the fit result onto the m(a,| ) (left) and m^-^^ (right) 



B{B 



■ peu) 



(1.08 ± 0.28(stat.) ± 0.16(sys.)x)10-'^, 
(2.57 ± 0.52(stat.) ± 0.59(sys.)) x 10"^. 



(17) 



We also set 90% C.L. upper hmits on the B rjlv , B r/'iu , B a^lv , and 5° 
aQiv branching fractions: 



B{B- 
B{B° 

9 Conclusions 



B{B- 
B{B- 
aliv )Bia\ 



< 1.2 X 10"^ (90% C.L.), 

< 4.5 X 10-^ (90% C.L.), 
ry7r°) < 5.3 x 10"^ (90% C.L.) 



)B{al r/7r+) < 3.3 x 10^^ (90% C.L.) 



(18) 



Preliminary results on charmless semileptonic B decays have been obtained on a sample of 88 
million T(4S') BB decays collected by the BABAR experiment at the PEP-II asymmetric-energy 
B factory at SLAC, in events in which one B meson decaying to a hadronic final state is fully 
reconstructed and the semileptonic decay of the second B meson is identified by the detection of a 
charged lepton. 

We have explored several new approaches to the study of -B ^ ^u^v decays, most of which are 
statistically limited and some of which are already competitive with the best existing measurements. 
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From the measurement of the spectrum of the invariant mass of the hadrons X„ (Mx) we derive 
the branching fraction 

B{B XJu) = (2.53 ± 0.29(stat.) ± 0.26(sys.)ia4i (theo.)) x IQ-^. (19) 

Prom the two-dimensional distribution of Mx and q^, the squared invariant mass of the two leptons, 
we derive the partial branching fraction for Mx < 1.7GeV/c^ , (f' > 8GeV^/c^ to be 

AB(5 ^ XJv) = (0.88 ± 0.14(stat.) ± 0.13(sys.) ± 0.02(theo.)) x 10"^ (20) 

Prom these two measurements, utilizing the calculation of the fraction of events in the selected 
phase space region from DPN [5] and BLL [6] in the Mx and Mx-q^ analysis respectively, we can 
extract 

\Vub\Mx = (4.77 ± 0.28(stat.) ±0.28(sys.)tg:|^(theo.)) X 10~^ and (21) 
\Vub\Mx-q^ = (^-92 ± 0.39(stat.) ± 0.36(sys.) ± 0.46(theo.)) x 10"^, 

respectively. These results are competitive with the current world average "^iVub] = (4.57 ± 0.61) x 
10-3. 

The theoretical error is dominated by the uncertainty on the shape function parameters. In the 
process of writing this paper we were informed of an improved determination of these parameters 
from the photon energy spectrum measured by Belle in 6 — 57 decays (see Sec. 2.1). Utilizing 
these new constraints on the shape function parameters we get 

{VublMx = (5.22±0.30(stat.) ±0.31(sys.)l[!:|(theo.)) X 10~3 (22) 
iVuhlMx-q' = (4-98 ± 0.40(stat.) ± 0.39(sys.) ± 0.47(theo.)) x lO'^, 

to be compared with the measured values of \Vub\ in Eq. 21. These results show that the deter- 
mination of the shape function parameters is critical to these measurements. Furthermore the 
two approaches deal differently with the shape function parameterizations: in particular the BLL 
approach does not correct the acceptance for shape function effects but assigns a large error. The 
corresponding result is therefore more stable but returns a worse error than the DFN approach in 
case of accurate determinations of the shape function parameters. 

On the other side, the dependence of the \Vub\ measurement as a function of the requirement 
on q"^ (Fig. 8) shows that both models reproduce the signal distribution well. 

The reconstructed Mx spectrum is also used to unfold the Mx distribution for B — > X^fi/ 
events and to measure its first and second moments (see Table 1). With a higher statistics data 
sample, these measurements can be used to constrain the shape function parameters. 

From the same data sample, several exclusive charmless semileptonic B decays are identified and 
their branching fractions measured (see Sec. 8.4). Imposing isospin and quark-model relationships 
we derive: 

B{B neu) = (1.08 ± 0.28(stat.) ± 0.16(sys.)) x 10"^ , 
B{B peu) = (2.57 ± 0.52(stat.) ± 0.59(sys.)) x 10"^ . 

Although these measurements are still statistically limited, the signal samples are very clean and 
are selected with very loose criteria which will allow the measurement of form factors with some- 
what larger data samples. We are grateful for the extraordinary contributions of our PEP-II 

'^See the HFAG page http://www.slac.stanford.edu/xorg/hfag/semi/winter04/winter04.shtml . 
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Department of Energy and National Science Foundation, the Natural Sciences and Engineering Re- 
search Council (Canada), Institute of High Energy Physics (China), the Commissariat a I'Energie 
Atomique and Institut National de Physique Nucleaire et de Physique des Particules (Prance), the 
Bundesministerium fiir Bildung und Forschung and Deutsche Forschungsgcmcinschaft (Germany), 
the Istituto Nazionale di Fisica Nucleare (Italy) , the Foundation for Fundamental Research on Mat- 
ter (The Netherlands) , the Research Council of Norway, the Ministry of Science and Technology of 
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dom). Individuals have received support from CONACyT (Mexico), the A. P. Sloan Foundation, 
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A The Unfolded Spectrum and the Covariance Matrix 



Table 7: Fraction of B ^ Xy^iu events measured in a given bin i of Fig. 6 and corresponding 
breakdown of errors. All numbers are in percent. 



i 


^unf 


a' 


•^stat ® "^MCstat 


^dct 


^n\v ® "^theo 


'^bks ® ^breco 


1 


2.890 


1.655 


1.206 


0.589 


0.322 


0.423 


2 


1.067 


1.294 


0.895 


0.415 


0.685 


0.266 


3 


5.355 


4.039 


3.325 


0.906 


0.952 


1.026 


4 


9.840 


4.934 


3.532 


1.388 


1.999 


1.537 


5 


28.024 


8.330 


5.302 


2.807 


3.111 


4.183 


6 


25.830 


3.657 


1.850 


1.223 


1.169 


2.377 


7 


15.880 


5.996 


3.964 


1.806 


2.742 


1.872 


8 


7.189 


5.742 


3.935 


1.940 


1.161 


2.998 


9 


2.742 


4.213 


2.534 


1.276 


1.939 


2.209 


10 


0.880 


1.822 


1.194 


0.627 


0.290 


1.097 


11 


0.261 


0.737 


0.478 


0.260 


0.104 


0.451 


12 


0.038 


0.149 


0.088 


0.064 


0.036 


0.084 


13 


0.002 


0.018 


0.006 


0.011 


0.008 


0.006 


14 


0.001 


0.005 


0.002 


0.003 


0.002 


0.002 


15 


0.000 


0.001 


0.000 


0.000 


0.000 


0.000 


16 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 
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Table 8: Full covariance matrix on the fraction of B ^ X^iv in the unfolded Mxspectrum (as in Table 7)). All numbers are in 10 
A version with a higher number of digits can be obtained from the authors. 



2.739 


1.021 


4.254 


5.685 


8.257 


-1.605 


-7.474 


-6.267 


-4.099 


-1.734 


-0.647 


-0.117 


-0.010 


-0.003 


-0.000 


-0.000 


1.021 


1.674 


2.975 


3.853 


4.058 


-0.992 


-4.398 


-4.012 


-2.465 


-1.153 


-0.460 


-0.091 


-0.009 


-0.003 


-0.000 


-0.000 


4.254 


2.975 


16.314 


13.999 


17.275 


-5.302 


-17.648 


-15.958 


-9.264 


-4.448 


-1.802 


-0.351 


-0.033 


-0.009 


-0.002 


-0.001 


5.685 


3.853 


13.999 


24.343 


34.115 


-1.585 


-28.030 


-24.898 


-17.030 


-7.222 


-2.692 


-0.484 


-0.039 


-0.012 


-0.002 


-0.001 


8.257 


4.058 


17.275 


34.115 


69.389 


11.837 


-45.085 


-45.370 


-33.773 


-14.300 


-5.325 


-0.965 


-0.082 


-0.025 


-0.004 


-0.002 


-1.605 


-0.992 


-5.302 


-1.585 


11.837 


13.373 


-0.222 


-5.212 


-6.439 


-2.675 


-0.980 


-0.176 


-0.015 


-0.005 


-0.001 


-0.000 


-7.474 


-4.398 


-17.648 


-28.030 


-45.085 


-0.222 


35.949 


31.214 


22.623 


9.146 


3.289 


0.574 


0.045 


0.014 


0.002 


0.001 


-6.267 


-4.012 


-15.958 


-24.898 


-45.370 


-5.212 


31.214 


32.971 


22.363 


10.272 


4.042 


0.763 


0.067 


0.020 


0.003 


0.002 


-4.099 


-2.465 


-9.264 


-17.030 


-33.773 


-6.439 


22.623 


22.363 


17.751 


7.207 


2.603 


0.466 


0.041 


0.012 


0.002 


0.001 


-1.734 


-1.153 


-4.448 


-7.222 


-14.300 


-2.675 


9.146 


10.272 


7.207 


3.319 


1.309 


0.249 


0.022 


0.007 


0.001 


0.001 


-0.647 


-0.460 


-1.802 


-2.692 


-5.325 


-0.980 


3.289 


4.042 


2.603 


1.309 


0.543 


0.107 


0.010 


0.003 


0.001 


0.000 


-0.117 


-0.091 


-0.351 


-0.484 


-0.965 


-0.176 


0.574 


0.763 


0.466 


0.249 


0.107 


0.022 


0.002 


0.001 


0.000 


0.000 


-0.010 


-0.009 


-0.033 


-0.039 


-0.082 


-0.015 


0.045 


0.067 


0.041 


0.022 


0.010 


0.002 


0.000 


0.000 


0.000 


0.000 


-0.003 


-0.003 


-0.009 


-0.012 


-0.025 


-0.005 


0.014 


0.020 


0.012 


0.007 


0.003 


0.001 


0.000 


0.000 


0.000 


0.000 


-0.000 


-0.000 


-0.002 


-0.002 


-0.004 


-0.001 


0.002 


0.003 


0.002 


0.001 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


-0.000 


-0.000 


-0.001 


-0.001 


-0.002 


-0.000 


0.001 


0.002 


0.001 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 
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